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SUMMARY 


A finite trigonometric series is fitted by harmonic 
analysis as an approximation function to the \!/-functIon 
of the Theodorsen arbitrary-airfoil potential theory. 

By harmonic synthesis the corresponding conjugate trigo- 
nometric series is used as an approximation to the 
e -function. a set of coefficients of particularly simple 
form is obtained algebraically for direct calculation of 
the e-values from the corresponding set of ^-values. The 
formula for e is 


n 



where the summation is for odd values of k only and 


v|/ 


k 




INTRODUCTION 


In the determination of the flow about an arbitrary 
airfoil (references 1 and 2) the problem arises of trans- 
forming a curve, nearly circular, into a circle. This 
transformation, a basic problem in conformal mapping, is 
further reduced to the determination of the following 
two conjugate Fourier series: 
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i ~ a o + / (%! cos mc P + s ^ n mc p) 


m=l 




e = (s m sin mcp - b m cos mcp) 


m=l 


( 1 ) 


(See references 1 and 2 for significance of notation.) 
fuo following integral relations are equivalent to series 
relations ( 1 ): 


P2tt 


'\i ( w ) = 


2 rr 


P2rr 


( cp T )dco ' + 


:TT 


<->0 


, V Qi » — Cp 

e(cp') co u — dcp' 


uo 


> ( 2 ) 


e ( r o ; 


1 _ 

* 2 rr 


' 2 tt 


\l/(cp’ ) cot ^ — 2 _i 2 dcpi 


^0 


It is convenient to introduce a new variable 
s = cp 1 - cp in relations (2). Because of the cyclical 
nature of these functions, the limits of integration may 
be written -it to rr. When the integral is broken into 
two parts, -it to 0 and 0 to rr, and -s is substituted 
for s In the first part, the following relations are 
obtained ; 


p2rr 

'l(cp) = ~ | \i/(cp)dcp + — 

err i d.,i 

°0 


n tt 


O 0 


e( r P + s) - e(cp - a)] cot y dr 


O TT‘ 


> (3) 


e(cp) = 


11 U 


— s 

v' ( cp — s) - \|/( cp + s) cot 5 - ds 


0 


J 


Thus, by use of relations (1), (2), or (3), e may 
be determined If \!/ is known or \j ; may be determined 
if e is known. 
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In the airfoil problem \|/ is specified as a 
function 1 of <p by means of a curve and e is to be 
determined. In theory the Fourier coefficients may be 
determined in relations (1) but in practice, because of 
the unknown analytic nature of the curve, it is neces- 
sary to resort to some type of numerical approximation. 

In references 1 and 2 an approximate method of han- 
dling the integrals of relations (2) is presented. In 
reference 3 a refinement of this method is given for the 
same integrals. An alternative procedure is to approxi- 
mate relations (1) by a finite trigonometric series and 
then to determine the coefficients by harmonic analysis. 
A development of this method is now given. 


HARMONIC ANALYSIS 


The ii-f unction is t o be approximated by a finite 
trigonometric series given by 


\i/ (cp) = Aq + A-j cos cp + . . . + A tl _ 1 cos (n - 1 )cp + A n cos ncp 

+ sin cp + ... + B n _]_ sin (n - l)cp 
n~l 

= A 0 + ) _(An cos mcp + Bpj sin mCp J + ^ cos ncp 
m=l 


If \i/ is specified at 


the range 
( 2n - 1 ) rr 
n 


0 ^ cp 1 2 it 
- then 


2n equally spaced intervals 

, . n it 2n 

- that is, 0, — , — , . . . 

n n 


in 


2n-I 


An = 


2n 


I 




r=0 


A In practice, \|/ is given as a function of 
0 = cp - e and therefore \| / is taken as a function of cp 
as a first approximation. An iteration process is neces- 
sary to determine both i s and e correctly as func- 
tions of cp. 
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where \[g 


= Hir) 


2n-l 


'V = n > t 


cos m 


rrr 

n 


r=0 

n--l 


Rfp 


\!/ s in m 


PIT 

~n" 


2n-l 


= 


dn 


(-D r X 


r=0 


B n = 0 


Now 


n-1 


€ ( cp ) = \ ^A m sin mg - cos mcpj + A n sin ncp 


i=l 


= /sin rnCo\~ 


r r*rr \ r> tt 

„ , .sin mcp) ^ cos m ; — > cos imp \ i]/„ s in re — 

i 1 / \ / 1 n X 1 n 

m ^ r r 


+ ~ sin ng^> (-l) r ii 

2 n / r 

r 


1 \ \ 

n /_ /_ 

m ' r 


X sin m 


(<P - ~) + sin ng ) (-lT % 


When the order of summation is interchanged. 
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i 5~ n ~^- , n ~-^- . 2n^_l 

e(<p) 35 n /> ^r/ sin m(cp - ~~) + — • sin ncp (-1)" 


r=0 m=l 


r=0 


Now if e is evaluated at the same points cp at which 
the values of \j/ are given, that is, at the points 

cp = E-Ul, the variable becomes 
n 

rTr - 1 - r )rr _ krr 

n _ ”n 


and the last term becomes zero. After this substitution 
is performed. 


\lr 

v r 



f v ’ tt krr\ 
v n n / 


' t '( cp + m) = 


and 


e (cp) 


2n- l n-1 
k=0 m=l 


m 


kTT 

n 


The summation with respect to k is also from 0 to 

2n - 1 because of the periodicity. But when k is odd. 


n-1 


s m m : 


ktr 


m=l 


n 


^ kir 
cot — 
2n 


and when k is even, 

n-1 


\ . k-rr 

y sin m — 


m=l 


n 


For odd values of k, therefore 

2n-l 

e (cp) = - - y \j; k cot 
k=l~ 


kTT 

2n 


0 
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/• 


t> 


or 


e ( cp) 




cot 


kTT 

2n 


Finally, then. 


e(cp) 


n 

_ V~ 

~ /_ 
k=l 


c k(*- k 



(i+) 


where 



^(c? 


+ 


krr^ 
n ' 


and for odd values of k 


C 


k 



krr 

2n 


and for even values of k 


C 


k 


0 


Equation (4) thus gives the same result as Is 
obtained by harmonic analysis and synthesis. Comparison 
with equations (J) indicates that equation (lq) may also 
be interpreted as the evaluation of this integral by the 
ordinary rectangular summation formula using intervals 
of width 2-rr/n and using the value of the integrand at 


the midpoint of each interval, that is, 
where k is odd. 


at 


kn 

n 


practical observations 


Squat ion (h) uses only one -half the available 
information. It is evident that all the points may be 
used because all the given n points may be considered as 
alternate (odd) points of a system of 2n points. The 
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vnlue of € so computed is, of course, to be plotted at 
cp-points midway between the given \|/ -points . The n values 
of ^ i / therefore give values of c corresponding to an 
approximation function consisting of a trigonometric 
series of n - 1 harmonics. Values of the coefficient C 
for n = 10, 20, I 4 .O, and 80 are given in table I. For 
smooth curves the present method for n = 20 is more 
accurate than the IpO-point method of reference 3 and 
requires only one -half as much computational virork. 

How to handle small irregularities or bumps in the 
curve is of interest. One procedure is to_fair through 
the bump and_to designate the faired curve \f/ . _The devi- 
ation from \\t is a Ad-curve. The conjugate e may be 
determined in the usual manner and a conjugate Ae may 
be determined by use of a very small interval, say, 
n - 200. The desired e-values are given by the sum of 
£ and Ae. This method cannot be justified on strict 
mathematical grounds but is probably more than adequate 
for engineering purposes. 


Langley Memorial Aeronautical Laboratory 
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TABLE I.- VALUES OP C k FOR USE WITH EQUATION ( 4 ) 


* 

c k 

1 

1 

n = 10 

n = 20 

3 

1! 

4=" 

O 

2 

II 

rr 

b 

1 

O.63138 

0.63531 

0.63629 

0.63654 

3 

.19626 

.20827 

.21122 

.21198 

s 

. 10000 

.12071 

.12568 

.12691 

ry 

i < 

1 9 

• 05095 
.01564 

.08159 

•05854 

.08364 

.06777 

.09037 

.06999 

ll 

13 


. 04270 
.03064 

Ml 

.05697 

.04790 

15 


. 02071 

.03742 

.04121 

17 


.01200 

.03171 

.03609 

19 


.00394 

.02704 

.03194 

21 



.02311 

.02858 

23 



.01971 

.02577 

25 



.01670 

.02339 

1 27 



.01400 

.02153 

29 



.01153 

.01955 

31 



.00922 

.01.794 

33 



.00705 

.01651 

35 



.00497 

.01523 

37 



.00296 

.01407 

39 



.00098 

.01500 

41 




.01202 

43 

H 

47 




.01111 

.01026 

.00946 

49 




.OO67I 

51 




. 00800 

53 




.00733 

53 




.00668 

57 




.00606 

CO 
j > 




.00547 

61 




.00489 

65 




.00433 

65 




.00579 

67 




.00326 

C 0 
L ' y 




.00274 

71 




.00223 

73 




.00173 

75 




. 00125 

77 




.00074 

79 ] 

.—I... 


}.... . — — — . — _ 

.00025 
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